Existence and stability of static shells for the 
Vlasov-Poisson system with a fixed central point 

mass 



Achim Schulze 
Mathematisches Institut der Universitat Bayreuth 
D 95440 Bayreuth, Germany 



March 13, 2008 



Abstract 



We consider the Vlasov-Poisson system with spherical symmetry and 
an exterior potential which is induced by a point mass in the center. This 
system can be used as a simple model for a newtonian galaxy surrounding 
a black hole. For this system, we establish a global existence result for 
classical solutions with shell-like initial data, i.e. the support of the density 
is bounded away from the point mass singularity. We also prove existence 
and stability of stationary solutions which describe static shells, where we 
use a variational approach which was established by Y. Guo and G. Rein. 

1 Introduction 

In stellar dynamics, the evolution of a large ensemble of particles (e.g. stars) 
which interact only by their self-consistent, self-generated gravitational field, is 
described by the Vlasov-Poisson system 



Here / = f(t,x,v) > is the phase-space density of the particles, where tEW. 
denotes time, and x,v S 1R 3 denote position and velocity. U — U(t,x) is the grav- 
itational potential of the ensemble, and p = p(t,x) is its spatial density. 

We want to investigate this system under the influence of a fixed point mass. 
If we assume that a point mass M c is fixed in the origin and acts like an external 
force on the particles, the Vlasov equation reads 



d t f + v-V x f-V x U-V v f = {), 
AC/ = 47rp, lim U(t,x) = 0, 



(1.1) 
(1.2) 




(1.3) 



d t f+v-V x f-(V x U-V xT ^ 



)-v /=o. 



(1.4) 



If we write U c s ■ = U — M c /\x\, the Poisson equation becomes 



AU eS = 4ir(p + M c 5), (1.5) 

where S denotes the (5-distribution. Here, we want examine to the existence 
and stability of steady states of the system Ql.2p - p.4p . One easily verifies 
that / = f(t,x,v) is a solution of (|1.4p . iff it is constant along solutions of the 
characteristic system 




V x (U(s,X)-m ' ^ 



where (X,V) = (X,V)(s):=(X,V)(s,t,x,v) with (X,V)(t,t,x,v) = (x,v) for an 
initial value (x,v)gM. 6 and s,teK. Thus for the construction of stationary 
solutions, a natural idea is to find conserved quantities of (|1.6p - now with time- 
independent U. One immediate expression for such a quantity is the particle 
energy 

and if we make some additional symmetry assumptions on the potential U, we 
can find other conserved terms such as the angular momentum. 

In this paper, we are interested in stationary solutions of the form 

fo(x,v) = (E -E) k + (L-L ) l + , (1.7) 

where 0<k<l, (•)+ denotes the positive part and E < 0,L > are constants. 
E is the particle energy as above and 

L=\xxv\ 2 = \x\ 2 \v\ 2 -(x-v) 2 (1.8) 

denotes the modulus of angular momentum squared which is conserved along 
characteristics, if U is spherically symmetric. 

If we want to construct the stationary solution (fo,Uo) explicitely from the 
ansatz (|1.7p . we still have to solve the Poisson equation (|1.2p to get a self- 
consistent potential IV The existence of stationary solutions with parame- 
ter range k >- 1 , 1>-1, k + l + l/2>0, k<3l + 7/2 was established in [TT] for 
M c = 0. 

Without the exterior potential, the existence and stability of stationary so- 
lutions of the form (|1.7p was done in [17], where the parameter range I > — 1 and 
< /c < Z + 3/2 was covered. 

As mentioned above, for our ansatz (|1.7p we require that the corresponding 
potential U is spherically symmetric and therefore the stationary solutions (|1.7p 
also have to be spherically symmetric, i.e., 

f(x,v) = f(Ax,Av) VAeO(3), (1.9) 

where 0(3) is the group of orthogonal 3x3 matrices. For Lq > the support of 
the induced spatial density p(x)=p(\x\) is contained in some interval [i?i,i?2] 
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with i?i > and the steady state describes a shell. This can be seen as follows. 
If we introduce the new coordinates r:—\x\,w:—x-v/r and L as in (|1.8p . we 
can calculate the spatial density of f as 

Pfo( x )= / fo(x,v)dv 

= ^JJ o X { E o-^ 2 + ^)-Uo(r) + ^j (L-L )' + dwdL (1.10) 
= C(k,l)r 2l ^E a -U Q (r) + ^-^) , (1.11) 

where 

C(k,l) = 2 l+ V 2 iT f 1 ^=ds f ' s l+x ' 2 {l-s) k ds. 
Jo VI -s Jo 

For small r the expression in the bracket of p. lip becomes negative and this 
implies suppp/ C [i2i,oo[, for some i?i>0. On the other hand, because of 
U' Q {r) — J Q r s 2 p(s)ds/r 2 >0 for r>0, the function —Uq (r) +M c /r is decreasing 
and with Bo < we conclude that pf (r) = for large r. 

These shells together with the exterior potential induced by a point mass 
can be used as a simple model for a galaxy which encloses a black hole in the 
center. 

The ansatz (11.71) also leads to steady states and shells of the Vlasov-Einstein 
system, the general relativistic counterpart of the newtonian Vlasov-Poisson 
system, and they provide an access to study stability and critical phenomena 
numerically, cf. [T]. 

We examine the shells in the newtonian framework and to investigate their 
stability we will firstly prove a global existence result for the system (jl.2l) 
(|1.4p for initial data, which vanishes in a neighbourhood of the singularity r = Q. 
The corresponding solution then exists for all time, and will always vanish, if 
a; is in a ball around the singularity, which is determined by the initial datum. 
We mention that, without the exterior potential, the existence problem for the 
Vlasov-Poisson system is well understood, see for example O [10l [16] for global 
existence of classical solutions. However, in our situation the exterior potential 
becomes unbounded in r = and we have to ensure that the particles stay away 
from the singularity. 

To show existence and stability of the shells, we use a similar approach as 
in [5] , where existence and stability of the above steady states was shown in the 
case Lo =0 without the exterior potential. The main idea is to use an Energy- 
Casimir functional as a Lyapunov function with the help of variational methods. 
Concerning this approach for stability issues for the Vlasov-Poisson system we 
also want to list [31 SI O El EH 03] here. We briefly sketch the basic concept: 
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The Vlasov-Poisson system is conservative, i.e., the total energy 
H(f):=E khl (f) + E pot (f) 

■=\ J \v\ 2 f(x,v)dvdx-±- J (w Uf {x)\ 2 + ^^p f {x)\dx (1.12) 

of a state / is conserved along solutions and hence is a natural candidate for 
a Lyapunov function in a stability analysis; Uf denotes the potential induced 
by /, note also the interaction term J pf M c /\x\dx induced by the fixed central 
point mass. However, the energy does not have critical points, but for any 
reasonable function <!> the so-called Casimir functional 




<&(f(x,v))dvdx 



is conserved as well. Now one tries to minimize the energy-Casimir functional 

H C :=H + C 

in the class of allowed perturbations Tm, which consists of positive L 1 (R 6 )- 
functions with prescribed mass M, i.e. JJ f = M and with finite kinetic energy 
and a finite Casimir functional to ensure that Tic is well-defined. 

The aim is to prove that a minimizcr /o is a stationary solution of (|1.2[) - (II. 4|) 
and to deduce its stability. One of the difficulties is to show that the weak limit 
of a minimizing sequence in Tic, indeed is a minimizer. For this purpose, we 
will need that every function in the class of perturbations Tm vanishes on the 
set < L < L . 

We are only able to show stability against spherically symmetric perturba- 
tions, because our approach requires an L-dependence in the Casimir functional, 
more precisely, we define 

C(f):= I ${(L-L )+ l f(x,v))(L-L ) l + dvdx, (1.13) 

JR 3 

with 0<k<l as in (|1.7[) . $ convex, satisfying certain growth conditions, and 
this will be a conserved quantity for spherically symmetric / only. To simplify 
our presentation, we focus on the case 

which will lead to stationary solutions of the form (jTTTJ) . The Casimir functional 
then reads 

C(f):= [ f 1+1 ' k {x,v){L-L )- l,k dvdx. (1.14) 

JS. 3 

At one point we need a scaling argument, which gets complicated in the case of 
a translation in L in the Casimir-functional. Here we exploit the spherical sym- 
metry and use coordinates adapted to it: If / = f(x,v) is spherically symmetric, 
we have 

f(x,v) = f(r,w,L), 
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with r = \x\, w= and L as in (|1.10[) . see Section 1.4. Altogether, we want to 
minimize the energy-Casimir functional 

H c (f)=E kin (f)+E pot (f)+C(f), 

with i?kin,^ P ot from (|1.12p and C(f) as in (|1.14|) over the set 



T M ■= i f e L L (IT) I / > 0, / is spherically symmetric, J J f = M, 

£ kin (/)+C(/)<oc, f(x,v)=0 a.e. for 0<L<L o \. (1.15) 



See (|1 .9[) for the definition of spherical symmetry. 

This paper is organized as follows: In the next section, we prove a global 
existence result for the system (|1.2p ~ (jl.4p . Afterwards, we examine the vari- 
ational problem and we show that TLc is bounded from below in Section 1.3. 
Then we prove a scaling property and that the the infimum of TLc is negative 
in Section 1.4. In Sections 1.5 and 1.6 we show the existence of a minimizer 
and analyse its properties; it is a stationary solution, and it is nonlincarly stable 
against sperically symmetric perturbations. 



2 Global existence 

In order to prepare the stability analysis, we want to prove a global existence 
result for classical solutions to the system (|1.2|l - (|1.4|) . so that we know that 
solutions in a neighbourhood of the examined steady states exist. We want to 
prove the following theorem. 

Theorem 2.1 Consider the system i f l.S^) ~ [L~4\ ). Let be a spherically 

symmetric initial datum with f{x,v) — for L := \x x v\ 2 < Lq, where Lq>0 is 
given. Then the corresponding solution (f,U) exists for all time and there exists 
Rmin>0, such that f(t,x,v) = for \x\<R m i n , t>0, where the number R m i n 
only depends on M C ,L$ and f . 

Remark. Without the exterior potential, the global existence result was proved 
by J. Batt, cf. [2] and this was also the first global existence result for the 
Vlasov-Poisson system in space dimension three. In our proof given below, the 
main idea for proving the boundedness of the velocities, is due to E. Horst, cf. 
0. 

Proof. We fix an initial datum / g C], (R 6 ) with / > and we fix R, P with 

f(x,v) = for | x | >-R or \v\ >P. 

This implies f(x,v) = for \x\<y/L~o/P, since L = sm 2 (a)\v\ 2 \x\ 2 > Lq on the 
support of /, where a denotes the angle between x and v. 
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In the following, we will denote first partial derivatives with respect to x with 
V x and we will write d 2 for the second partial derivatives. We now consider the 
following iteration process to construct the classical solution. The Oth iterate is 
defined by 

f (t,z):=f(z), t>0,z€M 6 . 
If the nth iterate /„ is already defined, we define 

Pn-=Pf n -= [ fndv, U n :=U Pn :=-p n *- r r, U nte $ := U Pn - -yf 
Jr 3 \-\ N 

on [0,oo[xR 3 , and for L — \x x v\ 2 > Lo we denote by 

Z n (s,t,z) := (X n ,V n )(s,t,x,v) (2.1) 
the solution of the characteristic system 

X = V, V = -\7 x U nteS (s,X) (2.2) 

with Z n (t,t,z) = z, where we want to examine characteristics which start on 
the support of /. We claim that \X n (s,0,z)\ is bounded from below by a posi- 
tive constant for all s > 0, n € N, so that the right-hand side of the charcteristic 
system is well-defined for all time. Together with (|2.ip - (|2.2l) this leads to the 
definition 

f u „ f( z n(0,t,z)) for z=(x,v): \xxv\ 2 >L Q 

; n+1 {t,x,v).-^ Q 

for the (n + l)st iterate. Note that, due to sperical symmetry, L = \X xV\ 2 
is a conserved quantity of (|2.2p and that ||/ n (i)||i = HPnCOHi = ||/||i since the 
characteristic flow is measure preserving. We introduce some notations: 

Po(*):=A 

P n (t) :=sup||Vfe-i(s,0,z)| | zesupp/, 0<s<i, l<k<n\, neN, 

R^ in (t):=M{\x\, (i,d)£M 3 xM 3 | (x,v) €supp/„(s), 0<s<t} 
= inf j|X„_i(s,0,2)| | zesupp/, 0<s<tj, neN. 

Next we show 

P„(*)<Co, Kun(t)>Ci, teR+, neN, 

where Co >0 and C\ >0 only depend on M c , ||/||i, ||/||oo an( A L . We abbre- 
viate X n (s) :=X n (s,0,z) for some zGsupp/ fixed. Now fix t >0; we then have 
X n eC 2 ([0,t]) with 

X n (s) = - (m Pn (a, \X n (s)\)+M c ) ■ , (2.3) 
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where we used the spherical symmetry and defined 

m pS S i r ) : =47T / T 2 p n (s,T)d,T. 



To get suitable bounds for the right-hand side of equation (|2.3p . we firstly use 
[15], Lemma PI: 



m Pn( s ,\ X n(s)\) 
\Xn( S )\ 2 



= \V x U n (s,X(s))\ 
<3(2n)^p r M\\ /3 y n (t)\\li 3 

<4.3 1 /3 7r 4/3 ||/|| l/3 ||/|| 2/3p2 (<) 

=:«i*(t). 

Furthermore, since L n (s) := \X n (s) x V^(s)| 2 is constant in s, 



l*n(*)| 2 > 



> 



which implies 



We also have 



Altogether, 



\V n (sW\sm 2 (Z(X n ( S ),V n (s)))\ ~ P n 2 +1 (s)' 
X n {s) 



\Xn(sW 



<-rPn+i(t), 0<S<t. 



\Xn(s)\<- 



+ M C 



\X n (s)\ <C**min 



\Xn(sW 



where 



M r 



C* = C*(||/|| 1 ,||/|| oo ,M c ,L )=max( ||/|| 1+ M c ,-^ + k 



Now define £ n (s) := (X n (s))j for i€ {1,2,3} and 0<s<t. Then 

\L(s)\<g(m), 0<s<t, 

where 

5 (r):=C*min(l,P n 2 +1 (t)l, rel. 
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If Cn( s ) 7^0 on ]0,t[, we have 



in(t)-U0) < L(t)-in(0) £„(t) + £„(0) 



= 2 



<2 / \£ n {s)\g(tin(s))ds = 2 / 5 (r)dr 
Jo J£(.[o,t]) 

<2 / 3 (r)rfr = 8C*P„ +1 (t), 



and hence 



Ut)-L(0) <2V2^P^( 2 1 (t). 
If £,n{s) — for some s€]0,i[, we define 

S _:=inf{se]0,t[ | £„(s) = 0}, s+ :=sup{se]0,i[ | = 0} 

and the calculation made above implies 

in (*) - L (0) I < I in (t) - in (s+) \ + | £„(*- ) - (0) 

<4V2C^P n 1 / 2 1 (t). 
Since £„ = = (V^),, we conclude 

P„+i(t)<P + 4V6C T P„ 1 ( 2 1 (t) t>0, neN. 

and therefore 

P n (t)<C , neN, 
where Co only depends on H/HijU/HocLoj-Mc and we also have 



Now we can continue with the iterates and prove their convergence. We have 

f n eC\[0M^ 6 ), ||/n(t)||oo = ||/||oo, ll/nWI|l = ll/l|l, t>o, 



f n (t,x,v) = for \v\> P n (t) or \x\>R+ / P„(s)ds, or L < L 



and 



p„eC 1 ([0,oo[xM 3 ), 



47T , 



||Pn(t)||l = ||/||l, l|p™WI|co<y||/||coP„ 3 W, *>0, 

p„(i,a:) = for|a;|>P+ / P n (s)ds 



Wc define 

|| v X 

U n , eS (t)\\ 

min,oo • — 

sup I |V x f7ft ie fif(i,a;)| | < |x| <oo > 



and ||^J7 ni efr(i)||niin,oo is denned analoguously. 

Now choose Tq>0. We want to prove that there exists a constant C>0, 
which only depends on To,/,Lq and M c , such that 

1 1 V xPn min,oo <C, te[0,T o ], nGN. 

In the following, C > may change from line to line, but there is no dependence 
on t€ [0,T ] orneN. We have 

\V xPn+1 {t,x)\< [ vJf(Z n (Q,t,x,v))] dv<C\\V x Z n (Q,t,-)\\* x , 

J\v\<P n (t) 1 J 

where 

||V x Z n (0,i,-)||^:=sup{|V x Z„(0,t,*)| | z:Z(0,i,z)esu PP /}. 

Next, fixi,oel 3 , t e [0,T ] and write Z n (s) = (X n ,V n ) (s) := (X n ,V n ) (s,t,x,v), 
where we require that i?„(0) € supp/. Differentiating the characteristic system 
with respect to x, we get 

\V x X n ( S )\ < \V x V n ( S )\, \V x V n ( S )\ < ||^^n,eff(fl)||mta,oo|V x X n («)|. 

By integrating and noticing V x X n (t) — E , V x V n (t) = 0, we have 
\V x X n (s)\ + \V x V n (s)\ 

<1+ f (l + ||dX,cff(T)||mi„,oo) (\V x X n (T)\ + \V x V n (T)\)dT. 
J s 

Gronwall's lemma now implies 

|V x X„(s)| + |V x y n (s)|<exp f (l + ||^[/„,eff(T)|| min)0O )dr, 

Jo 

and thus 

||V x/ o n+ i(i)||oo<Cexp / \\d 2 x U n Mr)\\ min,oo 

dr. 

Jo 

A well known estimate for the Poisson equation then implies, cf. [15j. Lemma 
PI, 

||3X,cffWILin.oo<C(l+ f \\dlU n . cif (r)\Un.oodT). 

Jo 

By induction, 

||dX,eff(i)|U tl> c <Ce c * 



and thus ||9^Z7 ra ,eff(£)||min,oo < C. Now we show that the sequence (/„) converges 
to some function /, uniformly on [0,T ] xR 3 xl 3 . For neN and z€l 6 , 

\f n+1 {t } z)-f n {t,z)\<C\Z n {0,t,z)-Z n _ 1 {0,t,z)\. 

For < s < t, we have 



|X„(s)-X n _ 1 (s)|< / |K(r)-V„-i(r)|dT, 



|V n («)-V»_i(«)|< 



< 



| V, U n (t, X n ( T ) ) - V, C/ n _ ! (r, X„ (r) ) | 
+ |V ie I7„_i,aff(T,JC n (r))-V ie C/ n _i, e ff(r,Jr„_i(T))| 
'"V a! l7 n (T)-V a ,I7 n _i(T) 



-2||dX-i,eff(r) 



) |X„(r)-X n _ 1 (r) 



where we used the mean value theorem and the factor 2\X u (t) — X n -i(r)\ in 
the last line is an upper bound for the length of a curve which connects X n (r) 
with X„_i(t) (s <t <t) and avoids the critical area B^j^i Cq - note again that 

we have the inequality i?" in (i) > \/Xo/Co- 

Recalling || II min,oo < C, adding these estimates and applying Gron- 
wall's lemma, we obtain 

\Z n {8)-Z n -i{ 8 )\<C J WVvUnW-VvUn-iWW^dT 

<C / \\Pn{r)- Pn-l{T)\\ 2 J\\p n (T)- p n -i(T)\\\ 3 d,T 

Jo 

<C / \\Pn(T) ~ P„-i(t)\\ dr 



<C ||/n(r)-/„_i(r)|L , 



where the second inequality follows by splitting the expression 



|V*l7(a:)|< 



p( x ) 
\x-y\- 



dy< 



\x-y\<R 



Pix) 
\x-y\ 



:dy- 



\x-y\>R 



P( x ) 

\x-y\- 



dy, 



and then using Holder's inequality and an optimization in R > 0, cf. |15j . Lemma 
PI. 

Also note that the support of both p n {t) and f n {t) is bounded, uniformly in n 
and t£ [0,To]. Altogether, we have 

\\fn+i(t)-fn(t)\\<C* f |I/„(t)-/„-i(t)|L dr, 
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and by induction, 



l/wW-Zn-iWlL^C-V^C-p, neN , 0<t<T o . 

I MOO n \ n \ 



This implies that the sequence (/„) is uniformly Cauchy and converges uniformly 
on [0,Tq] x K 6 to some function / G C([0,T o ] x M 6 ), which has the following prop- 
erty: 

f(t,x,v) = for \v\>C or \x\>R + C Q t. 

Furthermore, 

p n -*p~Pf, U n ->U:=Uf, (ti-kjo), 

uniformly on [0,To] x R 3 . Since To > was arbitrary, the proof is complete once 
we show that the limit function / has the regularity to be a solution to the 
Vlasov-Poisson system. With [T3], Lemma PI, we have 

1 1 V, U n (t ) - V x U m (t ) 1 1 oo < C 1 1 Pn (t ) ~ p m (t ) 1 1 H 3 1 1 Pn (t ) - Pm (t ) 1 1 

and 



\\%U n (t)-%U m (t)\\oo<C 



\\ Pn{t) - P,n(t)\\ oc 

+ d\\V x p n {t)-V x p m (t)\\ 00 +R- 3 \\p n (t)-p m (t)\\ 



for any 0<d<R. This implies that the sequences (V x Z7 n ) and (d%.U n ) are also 
uniform Cauchy sequences on [0,T ]xffi. 3 . Indeed, since all p n have compact 
support, uniformly in n, we can estimate 

IM*)-An(*)||l<C||/J n (t)-/V n (t)|| 00 <C||/ n (t)-/ in (*)|| 00 

which converges to zero. For the term with the derivatives of p n , we only know 
that 

1 1 V x Pn \ L) V2 Pm (t)\\oo<C 

with a not necessarily small constant C, but here we can choose d > in front 
of this term as small as we want. Hence we have 

V x U,d 2 x UeC([0,T }xR 3 ). 

Now we have for the characterstic flow Z, induced by the limiting field — W X U, 

Z= lim ^eC'^Tolx^rolxM 6 ), 

n — >oo 

and finally, 



/(*,*)= lim f(Z n (0,t,z))=f(Z(0 : t,z)), 

n — >oo 

so that /eC 1 ([0,T o ] xl 6 ) is a classical solution. □ 
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3 A lower bound on Tic 

We recall that we want to minimize 

U C (/) = £kin(/) + ^pot (/) +C(f), 

with E kin ,E pot from (II . and C(f) as in (|1.14j) over the set 
•7"M : = \f S L 1 (M 6 ) | />0, / is spherically symmetric, 



S kin (/)+C(/)<oo, f(x,v) = a.e. for 0<i<L 



Firstly, we want to establish a lower bound on Tic and we will need several 
estimates for pj and Uf induced by an element /eJF A/ . We will show that 
Epot(f) makes sense, that is, 

\7U f £L 2 (R 3 ) and \ ^-p s (x)dx < oo. 
Jr 3 \ x \ 



Lemma 3.1 Let n:=k + l + ^. Then there exists C>0, such that 
p) + "(xM~ 2l/n dx<C{C(f)+E^(f)), feT M . 
Proof. For any R> 0, we have 
Pf( x ) = / f{x,v)dv 



f(x,v)dv+ / f(x 7 v)dv 

\v\<R J\v\>R 

< I (L-L )f T f(x,v)(L-L Q f + 1 ^ T dv+-^ I \v\ 2 f(x,v)dv 



\v\<R 



R 2 



°{J\ |<fl (I ~ inl "' /r ' < / f 1+ Hx,v)(L-Lon l/k dv 
1 



R 2 



\v\ f(x,v)dv 



<C\x\^R^( f 1+ i(x,v)(L-L )- l/k dv 



R 2 



\v\ 2 f(x,v)dv. 



Optimization in R yields 



R: 



2 / \v\ 2 f(x,v)dv)\x\^r I / fi+i( x , v )(L-L )+ /K dv 



.-i/k 
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and thus 

<C\x\^h72 (E kin (f) + J f 1+ *(x,v)(L-L )- l/k dv^ ^ . 

Taking both sides of the inequality to the power 1 + — , dividing by r» and 
integrating with respect to x proves the assertion. □ 

From Lemma 13. II we see that a function / lying in Tm and its induced density 
Pf automatically are elements of certain Banach spaces which we now define: 

L k ' l (R 6 ) := {/:R 6 ^M measurable, spherically symmetric and 

|/| 1+ *(L-L )+ /fe efo ( /^ x 

equipped with the norm 

H/HmH // A 1 -* a -I.,,)-' Al.nl r 

and 

L n ' l (~R 3 ) :— |/?:R 3 — ->]R measurable, spherically symmetric and 
J \p\ 1+ i\x\' 2l / n dx<oo] 

with norm 



l|p||n,i:=(y \ P \ 1+ ±\x\- 2l / n dx 



Both spaces are reflexive Banach spaces. More precisely, / and pf are contained 
in the subsets Z/V(R 6 ) and L"''(R 3 ), respectively, which consist of the a.e.- 
nonnegative functions of these spaces. 

We now need some notations which clarify what E pot (f) and VC// means 
for / <^Tm- For spherically symmetric pGC^(R 3 ) Poisson's equation becomes 



l(r 2 C/'(r))Wp(r), 



where r:=|x| and we have £/sC 2 (R 3 ) and U'(r)—4:TrJ^s 2 p(s)ds/r 2 ) in par- 
ticular, WU(x) = U'(r)-. This motivates the following definitions. For /eJF^/, 
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i.e., p:=p f eL n > l (R 3 ), we define 

TOp(r) := / p(x)dx — 4:'K / s 2 p(s)ds. (3-1) 

J\x\<r JO 

U' p (r):=^f± (3.2) 

VW:=^f (3.3) 

/GO 

C/;( S )d S (3.4) 

and we will sometimes write [// or VJ7/ instead of C/ P/ = U p or WU Pf = V£/p, if 
p/ is induced by /. The definition (13. 3p implies 



/ |VC/ p/ (a;)| 2 da; = 47r / 



Now we can state the next lemma. 

Lemma 3.2 (a) Define the function ££C(K + ) by 



R* 1 for < R < 1 
i? 92 /or 1< R < oo ' 



where gi :=Z — fc+l/2>0 cmc? cfe :=4i + 5 — n > 0. Then there exists a con- 
stant C>0 such that for p(E L n ' l (M?) with J p(x)dx = M we have 

-E pot (p):= — J \VU p \ 2 dx + J -j^ P (x)dx 

1 r R ( m 2 (r) \ 1 

- 2 J 1 ~ J r^~ + 87TMc r/?(r) J dT + 2R + 2MMc ) 

< C((R) (1 + » ) + i (M 2 + 2Af M c ) , i? > 

where U p denotes the potential induced by p. 
(b ) For every R>0 the mapping 



T:L^(R 3 )3 P ^^\ [0tR] eL 2 ([0,R}) 



is compact, 
(c) For pi,p 2 eL n > l {m?)C\L 1 (M?) we have 



VU P1 ■ VU P2 dx = -Air I U pi pi dx 
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Proof. Obviously, we have m p (r) <M, and this shows the first estimate of (a). 
Now for pG L n ' k (R 3 ), we have 



rp(r)dr = 



L-k-l/2 2fc + 3 

v "+ 1 r "+ 1 p(r) dr 



< 



(r)dr 



<CR l - k+1 ^\\ P \\ n , 
< C #-*+V2 (1 + | H |i+*) 

where we used Holder's inequality in the second line. Furthermore, again by 
Holder's inequality, 



and thus 



K(r)|<Cr( 2 '+ 3 )/(" +1 >||p|| lM , r>0, 



fl !^ dr <C||p||2^+5-«)/(n + i) 



(3.5) 
(3.6) 



<C J R(«+ 5 -»)/("+ 1 )(l + ||p||'+«) ) 



which implies the estimate in (a). As to (b), by Q3.6P we already know that 
the operator T is bounded. To show the compactness of T, we use the Frechet- 
Kolmogorov criterion, cf. [18] . Theorem X.l. We take a bounded set K ei"' 1 
and to show the precompactness of TK, we redefine Tp := 7 -^rX[a.R] € £ 2 (R)- 
The crucial part is to show that 

IKTpJfc-TpHa-.O, fc->0, 

uniformly in p€K, where {Tp)h '■= (Tp)(- + h). For /i>0, we have 

m„(r + h) , , . m p (r) , . 



< 



h m 2 p (r + h) 



dr- 



v 

R-h 



h m 2 (r) f R m 2 (r) 



m 2 (r) 



■ dr 



R-h r 



R-h 



1 



\m p (r + h) — m p (r)\ dr 



(r + h) 21 

For the first four terms, one can use the estimate (|3.5p . Indeed, for example, 



m 2 Jr) 



R-h 



dr<C \\p\\ 2 nl A 4l + 4 ~ 2n y( n+1 Ur 



R-h 



= C\\p\\l l (R^ +5 ^^ n+1 ^-(R-h) 



(4/+5-n)/(n+ 
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and 



L m ? ir) —h-r dr = J h —{ 



rR-h m 2(„\ / u \ 2 



; 1 1 - '-"- h m=(r) f ft 

which converges to zero by Lebesgue's theorem. We have 



\m p (r + h)-m p (r)\<C\\p\\ n j (( r + fc)(2l+3)/(n+U _ r (2J+3)/(n+i)) ? 

and again by Lebesgue's theorem, also the last term converges to zero. Each 
term coneverges uniformly in p £ K and the case h < is completely analoguous. 
As to (c) , we firstly show the assertion for p\ , p 2 £ C°° nL n,! nL 1 . An integration 
by parts gives 



/ VU P1 ■ \7U P2 dx = 4tt [ U' {r)m P2 (r) dr 

J JTS.+ 



= 4irU pi {r)m P2 (r) 
= -47r / U Pl p 2 dx, 



r=0 



(4^) 2 / U Pl (r)r 2 p 2 (r)dr 



where the boundary term at infinity vanishes since |J7 Pl (r)| < ||/Oi||i/r and 
m p 2 ( r ) < IIP2II 1 and the boundary term at zero vanishes since m P2 (r) = 
0(r 2 ), r^>0. Now we consider approximating sequences {p\),{p2) C L n ' l C\ 
C°° PL 1 such that for i = 1,2 

pi^Pi in L"' 1 (i^oo), 

and || ||i < || pi ||i. Using the estimates of (a), we conclude that the above 
identity still holds for p t £ L"' Z nL 1 and the proof is complete. □ 



Lemma 3.3 There exists a constant C>0, such that 

Hc(f) > \ (E kin (f) + C(f))-C, f£T M 

in particular, 

h M ■= mi{Hc(f) I / G F M } > -00. (3.7) 
Proof. Using the previous two lemmas we have 

HcU) > E kin (f) +C(f) - C((R) (1 + \\p f ) - M2 + 2 fl MMc 
>(^ k i„(/) + C(/))(l-CC( J R))-CC( J R)- M '" ' ' ! U ' 



2R 

where C > is some constant which does not depend on R>0. The assertion 
follows by a suitable choice of R. □ 
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4 A scaling lemma 

In this section we show that Km is negative. We also examine the behaviour of 
Tic(f), if / is rescaled. 

Lemma 4.1 Define hu as in {3. ty - Then for M > we have — oo < Km < 0. 

Proof. As already mentioned in the introduction, we will use coordinates 
adapted to spherical symmetry. If f(x,v) — f(Ax,Av) V 0(3), we have 

f(x,v)=l(r,w,L), 

where r := \x\, w:= ^r, L := \x x v\ 2 and we will write again / instead of /. 

It is easy to check that, in the new coordinates, the energies and the Casimir 
functional read 

£W/) - 2n2 J R+ J & J M+ ( w2 + ^)f(r,w,L)dLdwdr, 

1 f m 2 t (r) f 
E V ot(f) = -o — ^j~-dr — 4wM e / rp f (r)dr, 
z Jr+ r Jm.+ 

C(/) = 4tt 2 / / / f 1+1 / k (r,w,L)(L-L )- l/k dLdwdr, 
Jr+ jr jr+ 



with R + := [0,oo[ and m/ =m Pf as in (|3.1| . 

Given any function /G-T-m, we define a rescaled and translated function 

f(r,w,L) = af (br, cw,b 2 c 2 L~ (b 2 c 2 - 1)L ) , (4.1) 

where a,b,c> 0. 

Then f(r,w,L) =0 a.e. if L<L , 



f(r,w,L)dLdwdr = a(bc) JJJ f(r,w,L)dLdwdr 
and if f GTm-, we have fGJ 7 ^ with M = a(bc)~ 3 M. Furthermore, 

E kin (f) = 2K 2 ab~ 3 c~ 5 [[[ (w 2 + L + {b2c l~ 1)Lo )f (r, w, L) drdwdL, (4.2) 



C(/) = a 1+ ^- 3 +fc- 3 +fc(/), (4.3) 



1 I .,2, «... fi "'/ (?jr) 



E P ot(f) = -^ I a^-Oc-o^^-dr-Airab- 2 ^ 3 I M c rp f (r)dr 



1 „2>-5 -6 / m /( r ) I ir ,;--2.. 



= --a 2 6^ 5 c^ b / -^—dr~4irM c ab~ 2 c' 3 rp f (r)dr. (4.4) 

To prove the lemma, we consider the case 6c < 1. Here we have 

E kin (f) < ab- 3 c- 5 E kin (f). (4.5) 
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Now we fix some / £ T\ with compact support and let 

a = M{bcf. 

Consequently, 

HcU) < a'+h- 3 +fc- 3+ fC(f) + ab-3 C ~ 5 E kin (f) 

--a 2 b- 5 c- 6 -J±-Ldr-AnM c ab- 2 c- z rp f (r)dr 
2 Jr+ f Jr+ 

<C 1 ai(bc)^+C 2 c- 2 ~C 3 b, 

where Ci,C 2 ,C 3 > depend on / and M. Since we want the last term to dom- 
inate as b~ >0, we let c = 6~' ) / 2 , so that 6c = 6 1 ~'5 for some ?yG]l,2[. For b small 
enough we have be < 1 and 

Hc(f) < db*- 1 -^ 21 ^^ + C 2 b v - C 3 b. 

Now fix r\ €]1,2[ such that (1 — §)(2Z + 3)/k > 1; such an r\ exists by the assump- 
tions on k and /. For 6>0 sufficiently small, the sum of the last three terms 
will be negative and the assertion follows. □ 
In the next section, we will use the rescaling formulas (|4.2p - (|4.4p to show that 
a function /o, constructed by the weak limit of a minimizing sequence actually 
is a minimizcr with mass M . 



5 Existence and properties of minimizers 

Theorem 5.1 Let M > 0, Lq > and let (fj) C Tm be a minimizing sequence of 
He- Then there is a minimizer /o and a subsequence (fj k ) such that Hc(fa) = 
Km and fj k — 1 fo weakly in L k ' 1 . For the induced potentials we have V{7j fc — > Vi/o 
strongly in L 2 (K 3 ). 

Proof. By Lemma [5751 i?kin(/j) +C(fj) is bounded and thus (fj) is bounded in 
L k ' 1 . Now there exists a weakly convergent subsequence, denoted by (fj) again: 

fj-^fo weakly in L k ' 1 . 

Clearly, fo>0 a.e. and fo(x,v) — a.e. for < L < Lq. By weak convergence, 

Ekin(fo) < limsup^kin(/j) < 00. (5.1) 

j— »°o 

By Lemma [37T1 (pj) = (pt.) is bounded in L n ' l (R 3 ). After choosing another 
subsequence, we conclude that 

Pj^Po weakly in (5.2) 

where we have the identity 

Pa=Ph ■= / fo(x,v)dv. 
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Indeed, assume we would have pf a > pq a.e. on the measurable set A :— :— 
{xGR 3 |i?i < |x| <i?2 with < Ri < i?2 < oo}, note that both p and p/ are 
spherically symmetric. Then for i?>0, by weak convergence we have 



0<7:= / (pf (x)-p (x))dx 



- lim 

i^°°J aJ\v\<r 



- lim / pj(x)dx, 



fj{x,v)dvdx + 



fo(x,v)dvdx- 

AJ\v\>R 



where we used the fact that \A S (L n ' 1 )* and xaxb r € ( y L k ' 1 )* . Now E^ n (fj) is 
bounded and this implies 

2 2 C 

f (x,v)dvdx< -^E kin (f ) <-^limsup£kia(/j) < 



'A J|w|>fl 

We conclude 



i? 2 



i? 2 



- lim 

j->°°JAJ\v\>R 



R 2 



fj(x,v)dvdx< 



R 2 ' 



2C 



(5.3) 



which is a contradiction. 

Next, from (|5.2|) together with Lemma I3T21 (a) (b), the strong convergence 

VUj -> VC/ strongly in L 2 (M 3 ) , 

follows, and we have 

Epot {fj ) — > Spot (/o) • 

Indeed, from Lemma 13.21 we have 

™2 



1 

-in 



\\7Uj-\7U Q \ 2 dx = 



< 



Pf, -Pfo 



dr 



R m 



ps 3 -pso , M 2 
r it 



:/+//. 



Now let e>0 be given. Choose R>0 large enough so that II<e/2. For j 
sufficiently large, the first term also will be smaller than e/2 because of the 
compactness of T, defined in Lemma 13.21 (b): The weak convergence pf. —^Pf 
implies the strong convergence m p/ _ P/o /r— >0 in L 2 ([0,i?]). 
Furthermore, we can estimate the interaction term as 



(p J (x)-p Q (x))dx 



< 



B R \X\ 



■(pj(x)-p {x))dx 



+ ■ 



2M 
~R 



Here the first term tends to zero, because of the weak convergence (|5.2p together 
with the fact that (]jj, - )l 2 (Bh) ^ (£™''(R 3 )) which we have shown in the proof 
of Lemma l3.2f a). The same argument as above then proves 



R 3 \ x \ 



■{pj(x)-p (x))dx 
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Next, we show that fo actually is a minimizer, in particular 
Ekin(fo) +C(fo) < oo. By weak covergence, we have 



C(/o) = ||/o|lT l; ^<lin 1 inf||/,|| K r /R <oo. 

3—>oc ' 

Together with (|5.ip and (|5.3|) this implies 



£ki„(/o)+C(/o)< lim (B k »(/ J -)+C(/ J -))<oo, 
note that the linx,—^ in the above inequality exists. Finally, 
«c(/o)=C(/o) + £kin (/o) + £ po t (/o) < lim (C (/, ) + B kin (/, ) + E pot (/,)) = /i M . 

J— >oo 

It remains to show that ||/q||i = M. By weak convergence, we have ||/o||i <M 
and we already know that ||/o||i>0, since fiM<0. Now assume that Mq:= 
\\fo\\i<M. We consider the rescaled function fo defined in (|4.I| in section 4 
and recall formulas (|4.2p ~ (|4.4p . Now define 



o:=1 ' c:= ImJ ' h:=c ■ 

This implies (6c)" 3 = M/M and thus ||/ || 1 = M. We have 
h M <H c (fo) 

<c£ kin (/ ) + c 3 - 2i / fe C(/ )--c 4 / -Z±J-dr-ctorM r Pf (r)dr, 

1 JR+ r JO 

where we used (|4.5[) . note that &c = c _1 < f . Since c > 1 and < fe < / we conclude 



h M <Hc(f )<cHc(fo)=[ j^J h M , (5.4) 
which is a contradiction. □ 



Theorem 5.2 Let /o6^m 6e a minimizer of Tic- Then there exists Eq<0 
such that 

f a (x,v) = ] -^ 1 (E -E)l(L-L ) l + (5.5) 

1 , , , M, . , 

E: =2 V + Uo{x) -J x J (5 ' 6) 

awe? J7o is the potential induced by fo. Moreover, fo is a steady state of the 
Vlasov-Poisson system il.S\) - (L4^ . 
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Proof. Let /o be a minimizer. We choose a suitable representative for /o and 
define for e > the set 

e < /o (a?, w) < i , i + e<i<io+^ 

Since / € i fc,z we have < \K e \ < oo for e sufficiently small. Now let g € L°°(R 6 ) 
be spherically symmetric with supp gCK e , and 

h ' =9 ~\E~\\jj 9dvdx ) ' XK '- 

Then for tGR small enough we have fo + rh>0 and /o+t/igJm, indeed, 
-Bkin(/o + T^)<oo and 



C(/o + r/ l ) = C(/ ) + r yy $'(/ )(L-L ) + ^ + O (r)<oo, 

where we recall that $(/) = / 1+1/fc . Now we have 
0<Hc(fo + Th)-H c (fo) = 

= r JJ (*\f )(L-L )- l/k + ±v 2 + U (x)~^hdvdx + o(T) 

= t JJ ($'(/ )(L-L );' /fe + E)hdvdx + o[r), 

where we used Lemma 13.21 (c) to calculate the potential energy term. Since —h 
is also an admissible function, this implies 

$'(/o)(L-X )+ /k +E]hdvdx = 0. 



Inserting the definition of h we get 

(&((L-Iv)+ l f )+E)-j±-ff (&((L-L )+ l f )+E) gdvdx = 0. 
Consequently, 

<f>'((L-L )+ l f )+E = E e a.e. on if £ , 

where 

Ee:= W]JJ (^'(( L - L o)+fo) + E)dvdx. 

Thus for e small, -E e will be a constant which we denote by Eq and we conclude 

$'((L-L Q )+ l f ) + E = E a.e. on {(x,v)\f (x,v)>0}. (5.7) 

Suppose now, there would exist a measurable set A C {(x,v)\fo(x,v) = 0,Lq < L} 
with 

E<E a.e. on A 
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and 0< \A\ <oo. We can also assume that A is spherically symmetric, i.e. \A 
is spherically symmteric. Next, define 

h:=XA ~J^~\(JJ XAdvdx^j -\K t 

with K e as above and small e>0. Then for r>0 sufficiently small we have 
/o +rh S Tm and again 

0<H c (fo + Th)-Hc(fo)=T JJ (<S>'((L-L Q )+ l h) + E)hdvdx + o{T). 

Plugging the definition of h into the above equation, we have 

0< II {&{{L-L )- + l f ) + E) X A-E Q JJxa 

{E-E )<0, 



a contradiction and thus E>Eq a.e. on {(x,v)\fo(x,v) = 0,Lo < L}. Together 
with (|5.7p this implies that /o is of the form given in the theorem. 

Since / is a function of the microscopic energy E defined by (|5.6p and L, it 
is constant along solutions of the characteristic system 

f X = V 

| V = -V X U (X)-^X 

and thus fo is a solution of the Vlasov equation, provided the potential Uo is 
sufficiently smooth. But one can indeed show that UoEC 2 (M. 3 ). This can be 
seen as follows. We firstly recall the formula for pf , if / is of form 



M j xfc+i+3/2 

Pair) := Pfo (r) = C(k,l)r" [ E - U (r) + - -£ (5.8) 



and we claim that £/o G £°°(;R + ) and thus the above equation implies pf S 
i^nL 00 . Indeed, for any i?>r, 

-U (r)=l R ^ds+r^ds 

J r S J R S 

<cj\^ k -^ + V\\p \\rud S +^ 

= CHpoIki (i?(- fc + i +V2)/(n+D _ r (-fc+i+l/2)/(n+l) J + ^£ 

and because of 0<fc<Z + l/2, the claim follows. Now / o/ €L 1 ni 00 implies 
Uq € C 1 and because of (|5.8[) also p n 6 C 1 . Together with £/o( r ) = Jq s2 Po(s)ds, 
the asserted regularity of C/q is proved. 
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By construction, we have 

so that {fo,po,Uo) is indeed a solution of the Vlasov-Poisson system. It remains 
to show that Eq<0. Recall the formula for po from (|5.8[) and the fact that 
|| /o ||i = M. If E > 0, we would have 

||/o||i = ||po||i>C(fc,0 J r r 2 <+ 2 (^j dr = C ^ r^ 1 / 2 dr = oo, 

where we have chosen i?o>0 sufficiently large so that Lo/r 2 <M c /2r, r> Rq. 
Consequently we conclude E < 0. □ 



6 Dynamical stability 

We investigate the nonlinear stability of fo . For / G Tm , 

Hc(f)-H c (fo) = d(f,f )-± J \WU f -WU fo \ 2 dx, (6.1) 

where 

d(f,fo):= II [(/ 1+1/fc -/o 1+lA ) (L-L Q )- l / k + (E-E Q )(f-f Q ) 



dvdx, 

where E is defined as in (|5.6[) . We have d(/,/o)>0, /gJFm with rf(/,/o) = 0, 
iff f = fo- Indeed, 

rf(/,/o)> ^ [$'((£- £ )+/o) + (£--Eo)] (f-fo)dvdx>0, 

which is due to the convexity of $, and on the support of fo the bracket vanishes. 
This fact allows us to use d(.,fo) to measure the distance to the stationary 
solution /o. 

Theorem 6.1 Assume that the minimizer fo is unique in Tm- Then for all 
e>0 there is 5>0 such that for any solution f(t) of the Vlasov-Poisson system 

with /(0)eC* c 1 (R 6 )nJ c 'A/, 

</(0),/o) + ^ I \\7U m -X7U fa \ 2 dx<5 



implies 



d(f(t)Jo) + ^ J \VU f(t) -VU h \ 2 dx<e, t>0. 



Proof. We observe that He is conserved along any solution f{t) of the Vlasov- 
Poisson system with /(0) € CH^^Tm ■ This follows from conservation of 
energy and the fact that both f(t) and L are conserved along characteristics. 
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Assume the theorem were false. Then there exists eo>0, tj>0, and fj(Q)£ 
C c 1 (R 6 )nJ r M such that 

d(M0)Jo) + ^ J |W /j(0) -V[/ /o | 2 dr<i 
d(fj(tj)Jo) + -^ J \VU fjitj) ~VU fo \ 2 dx>e . 



and 



From (16.11). we have 



lim Hc(fj{0)) = hM, 

J-tOO 

and because Hc(fj(t)) is conserved, 

lim H c (f j (t j ))= lim H c (f j (Q))=h M . 

J — >oo J — >oo 

Thus (fj(tj))cFM is a minimizing sequence of He and with Theorem 15.11 we 
have 

^ / |V[/ /3(tj) -W /o | 2 ^0, 

which implies 

d(/j(*i),/o)-0 

by (|6.1j) . a contradiction. □ 



Corollary 6.2 // m Theorem \6.1\ the assumption ||/(0)||k,i = ||/o||fc ; is added, 
then for any e > the parameter 5>0 can be chosen such that the stability 
estimate 

||/(t)-/o|kl<e> <>0 

holds. 

Proof. We repeat the proof of Theorem 16.11 except that in the contradiction 
assumption have 

\\m)-M\k,l + d(f j (t j )J ) + ± J \VU Mtj) -\7U fo \ 2 dx>e . 

From the minimizing sequence fj(tj) we can now extract a subsequence which 
converges weakly in L k ' 1 to /o. But due to our additional restriction we have 

||/i(*i)lki = ||/olki, JGN. 

Now the lower semicontinuity of the norm and the uniform convexity of L k ' l (R 6 ) 
imply fj(tj) — + fo strongly in L k ' 1 . Together with the rest of the proof of Theo- 
rem EUl the assertion follows. □ 
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Remarks. 



(a) The technical assumption / = a.e. for < L < Lq in the class of pertur- 
bations J-m, see (11.15[) . is needed for the scaling argument in Lemma [4~T1 
and it would be desirable to improve it to / = a.e. for < L < jL for 
some < 7< 1. 

(b) For M c = one can show existence and stability for steady states of form 
(fL~7)) for the parameter range I > -1, 0<fc<£ + 3/2. see [17]. ForM c >0, 
we had to restrict the parameter range to < k < I in order to guarantee 
that the scaling argument (|5.4|) works. 

(c) The uniqueness of the minimizcr /q subject to the fixed mass constraint 
can be shown by a scaling argument in the case Lq = and M c = 0. For 
Lq > 0, at least numerically the minimizer seems to be unique, but the 
scaling argument fails because of the translation in L. We mention that, 
for Theorem 16. 11 it would suffice if the minimizers of Tic were isolated. 

(d) We only obtain stability against spherically symmetric perturbations, be- 
cause the quantity L is conserved by the characteristic flow only for spher- 
ically symmetric solutions. Stability against asymmetric perturbations is 
an open problem and more delicate mathematical tools have to be invented 
to address this question. 
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